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We study weak itinerant ferromagnetism in one-dimensional Fermi systems using perturbation 
theory and bosonization. We find that longitudinal spin fluctuations propagate ballistically with 
velocity v m <C vp, where vp is the Fermi velocity. This leads to a large anomalous dimension in 
the spin-channel and strong algebraic singularities in the single-particle spectral function and in the 
transverse structure factor for momentum transfers q ~ 2A/vp, where 2 A is the exchange splitting. 
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Recently several authors presented conductance mea- 
surements in ultra low-disorder semiconductor quantum 
wires and suggested that an unusual feature in the range 
0.5— 0.7x 2e 2 /h of conductance can be explained in terms 
of spontaneous ferromagnetism |],§, §• At first sight this 
interpretation seems to contradict the Lieb-Mattis theo- 
rem B, which rules out magnetized ground states for 
electrons moving on a line, as well as for one-band lat- 
tice models in one dimension (Id) with nearest-neighbor 
hopping and interactions involving densities. However, 
there is no fundamental principle that forbids ferromag- 
netic ground states in quasi Id systems with finite width 
or one-band lattice models in Id with more general hop- 
pings. Indeed, numerical studies || show that the ground 
state of the Id Hubbard model with hopping between 
nearest and next-nearest neighbors can be ferromagnetic 
in a substantial range of densities and on-site interactions 
U. Clearly, the precise form of the energy dispersion 
e k plays an important role in stabilizing ferromagnetism 
H 7], ||. In principle, it should therefore be possible to 
design metallic systems with ferromagnetic ground states 
by properly adjusting the hopping integrals between the 
relevant orbitals. A promising class of Id materials where 
this might be achieved are certain types of organic poly- 
mers j^] , whose molecular structure can be designed in a 
controlled manner in the laboratory. Motivated by these 
new developments, in this work we shall use a combina- 
tion of perturbation theory and bosonization to derive 
some physical properties of itinerant ferromagnets in Id. 

Let us briefly consider this problem from a renormal- 
ization group (RG) point of view. The usual RG ap- 
proach to Id metals is based on the assumption that 
their long-wavelength and low-energy properties are de- 
termined by wavevectors k in the vicinity of the Fermi 
wavevectors ±kp. Given a general energy dispersion e k , 
it therefore seems reasonable to expand for k close to k F 



£fc F + v F (k - k F ) + 



(1) 



and similarly for k w — hp. By power counting, the Fermi 
velocity v F is a marginal coupling, while the inverse ef- 



fective mass 1/m* and the cubic parameter A are irrele- 
vant in the RG sense. In the field-theoretical formulation 
of the RG jl0| , these irrelevant couplings are simply ig- 
nored. However, as shown below, the cubic term in Eq. 
(|l|) is crucial to stabilize a ferromagnetic ground state 
in Id, so that a proper RG treatment of itinerant ferro- 
magnetism should include also the irrelevant couplings 
associated with band curvature effects. Therefore meth- 
ods which cannot properly handle these couplings, such 
as the field-theoretical RG |lCj] or bosonization, lose much 
of their power. Nevertheless, as shown below, in certain 
regimes bosonization is still useful to obtain nonpertur- 
bative results for correlation functions. 

We consider the following Hamiltonian describing in- 
teracting electrons on a Id lattice with length L, 



ka 



1 

21 



^2fijPi(-<i)Pj(<l) . ( 2 ) 



where ct and c ka are creation and annihilation opera- 
tors for electrons with momentum k and spin a. The 
labels i and j assume values in {n,m}, where n corre- 
sponds to the charge density p n {q) = Ylka ^la^k+qa' an d 

m denotes the spin density p m (q) = J2ka a ^ka^k+qa- To 
discuss spontaneous symmetry breaking we should start 
from a spin-rotationally invariant H, which constrains 
the bare to satisfy f nm = f mn = and precludes 
any momentum-dependence of f m = f mm . We also take 
/„ = f nn to be momentum-independent 111) . 

As a first step, we study the ferromagnetic instability 
within Hartree-Fock theory. Adding and subtracting the 
counterterm A (T (m) = f n n+af m m, where n — (p n (0)}/L 
is the density and m — (p m (0)}/L is the spin density, we 
may write H = Hq + Hi , with 

H - pN = ^J ka c ka - | [f n n 2 + f m m 2 ] , (3) 

ka 

and Hi = (2L)" 1 J2 q ,i fiSpi(-q)Spi(q). Here £ kt7 = e k - 
p, + A a (m) is the Hartree-Fock energy, Spi(q) — pi(q) — 
5q,o(pi(0)}, and N — /5„(0). In the ferromagnetic state 
the Fermi wavevectors k a and velocities v a are defined by 
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e h a — M + A CT (m) = and u^- = dek/dk\ k , while in the 
normal state Ck F — [i + A (T (0) = and vf = dek/dk\k F - 
Hence — ej, F + f n 5n = crA where A = — f m m and 
5n = n(m) — n(0). For convenience we keep the chemical 
potential fj, constant, so that the density n is a function 
of to. The two equations e^^ — ek F + f n 5n = oA, a — 
±1, together with the self-consistency conditions to = 
7T (fcf — fc|) and Sn — 7r _1 (fc| + fc| — 2kp) fix the four 
quantities fc-f, A^, <5n, and to. 

Throughout this work we shall assume m -C n (weak 
ferromagnetism) . The low-energy properties are then de- 
termined by wavevectors in the vicinity of the Fermi sur- 
face, as discussed in the classic work by Dzyaloshinskii 
and Kondratenko |fl2| . Hence we may expand around 
±kp. To leading order, it is sufficient to truncate the ex- 
pansion at the third order, see Eq. (|l|). Keeping in mind 
that 7tto <C &f and defining q m = A/vf we obtain 



2(1 + F ) 



2aAq\ 



(4) 



where A — j^{\2 — ), with Ai = \/{m*vp), A2 
X/vp, and Fq = 2f n /irvF- The Fermi velocities are 



A 2 

Va/v F = 1 + 0"Aig m + ^ I A 2 - - 1 

z V 1 + -TQ 



(5) 



Substituting Eq. ([I|) into m = 7r _1 (fc| — Aj), it is easy 
to see that, besides the solution m = 0, there is a non- 
trivial solution 7ttoq = 

[2(ii, - l)/(/^)] 1/2 , provided the 
radicand is positive. Here Iq = —2f m /TTVF is the di- 
mensionless Stoner parameter ||. To see whether the 
solution too is stable, we consider the energy change 
S£lo{m) = 17o(to) — f2o(0) due to a finite value of to, 
where SIq ( m ) = (Ho — /iiV) . We obtain 



<5£1 (to) = 



47T 



-Jo(/ - 1)(^™) 2 + j^™) 4 



(6) 



Obviously, a necessary condition for too to represent a 
minimum of Qo(m) is A> 0. In addition, the square root 
[2(I - l)/(J(f A)] 1 / 2 i s only real if either I < or I > 1. 
For consistency, we should also require that 7ttoo <C kp 
and that the band-structure is such that the higher order 
corrections in Eq. (^|) are small. For some special form 
of efc it should be possible to satisfy these conditions 
even for small negative Io provided hp A 3> 7r 2 |/o| -3 . 
Here we shall not further consider this case, but focus in- 
stead on the regime close to the Stoner threshold, where 
Io is slightly larger than unity. The distance from the 
critical point is then measured by the small parameter 
<5o = 2 (Jo — l)/^o- Interestingly, the numerical results of 
Rcf. U indeed show a critical Iq of order unity for not too 
large densities, which suggests that even in Id the Stoner 
criterion can be a reasonable estimate for the ferromag- 
netic instability. 



For simplicity we now set f n = — f m = /o > 0, cor- 
responding to a repulsive Hubbard on-site interaction 
p| . Note that close to the phase transition Iq = Fq = 
1 + O((5o). Let us first consider the density-density (xnn) 
and the longitudinal spin-spin (xmm) correlation func- 
tions. Within the random-phase approximation (RPA) 
we obtain 

XZ A (Q, ™) = [X? T + Xu - 4f X° n X Q u }/D , (7a) 
X ^(q, in) = [X° n + X° U + 4/oX? T X?J/£ , (7b) 
where D(q,iuj) = 1 - 4/ 2 X| T Xh, and 

/ - \ 1 f(£k+q/2,<T>) ~ f(£,k-q/2,a) , , 
Xaa'{Q>^) = ~T / 7 — • ( 8 ) 



L k £k+q/2M> 



S.k-q/2,, 



Here f(E) is the Fermi function. For small q and u) we 
may approximate 



xlM-> iuJ ) 



q 



7T (Vcrq) 2 + 



(9) 



For w > the dynamic structure factors Sf FA (q, w) 
ir~ 1 Imx^ PA {q, w + i0) can then be written as 



Sr A (q,u) = Z i \q\5(u;-v i \q\) 



(10) 



with Z n ■ 



-- [ttV 1 + F ] \ v n = vp^Jl + F , and Z m = 
v m — vfV$o~- Note that Sj iPA (q,uj) sat- 



F T L and 

The latter 
1 _ 



isfy the sum rules |lj 21im g ^ / °° % S? PA (q,u) = X i 
with the compressibility Xn = [ttvp(1 
the spin susceptibility Xm — 2/(7ru^(5o) 
is related to the Hartree-Fock energy (||) via x 
L^ 1 d 2 Qo(m) /dm 2 \ ma . We conclude that longitudinal 
spin fluctuations in Id can propagate ballistically, with 
velocity v m <C vf. In contrast, in 3d itinerant ferromag- 
nets the longitudinal spin mode can decay into particle- 
hole pairs and is therefore strongly Landau-damped || . 

Next, let us calculate the transverse spin-spin correla- 
tion function Yu(<?, ilo) within the ladder approximation 
shown in Fig. |l|, which yields 



X^f J (g,»w) = [x? 1 (?,iw)" 1 -2/o] 
For \q\ <C q m and \uj\ <C A we may expand 



Xu(9.*w)« 2A 



IU1 „ o 

2A-^ 



(11) 



(12) 



with the nonuniversal constant @B= jL[A 2 -A 2 ]. Note 
that B > A > 0. Using A = /o m o we obtain 



LAD (<z,H 



to 



bq 2 



(13) 



where b = 2AB is the spin wave stiffness. This im- 
plies a (5-function peak in the dynamic structure factor, 
S-\i(q,n) — too<5(w — bq 2 ), which exhausts the sum rule 
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FIG. 1: Ladder approximation for xtj.i see Eq. (|l l|). The solid 
arrows are the Hartree-Fock single particle Green functions 
with given spin projections and the wavy lines represent the 
bare interaction. 



FIG. 2: Diagrams contributing to Xtl m bosonization, see 
Eq. ([l5]). The thick wavy line represents the RPA interaction. 
To second order in the RPA interaction there are already 13 
diagrams contributing to xtl- 



la Sllil'^) — m o/bq 2 - The existence of well-defined 
transverse spin waves in the symmetry broken phase fol- 
lows from general hydrodynamic arguments |l3| |. How- 
ever, in Id it may well be that interactions lead to anoma- 
lous damping of spin waves and a breakdown of hydro- 
dynamics. This problem deserves further attention. 

Because the ferromagnetic instability is triggered by in- 
teractions with zero momentum transfer, we expect that 
at low energies the relevant interaction is dominated by 
forward scattering. Moreover, for m — it is known 
that repulsive backscattering interactions are marginally 
11T elevant Jk}. We assume that this remains true in the 
ferromagnetic state and expect that this assumption can 
be verified using RG methods. Note also that weak fer- 
romagnetism in 3d can be understood within the frame- 
work of Fermi liquid theory so that it is natural 
to expect that Luttinger liquid theory is the correspond- 
ing low-energy theory in Id, at least if the characteristic 
magnetic wavevector q m is small compared with hp. 

The leading long-distance behavior of correlation func- 
tions can then be obtained from a generalized Tomonaga- 
Luttinger model, where the energy dispersion is lin- 
earized around the Fermi points ±k a . Introducing a 
bandwidth cutoff A such that q m -C A <C &f and defin- 
ing field operators tp^(q) — \Lc a k^+q,cr, where a = ±1 
labels the Fermi points, the kinetic energy is represented 

by Eaa I- a If ^l^aHl^ail)- The interaction is for- 
mally identical with Eq. (^|), but with an implicit mo- 
mentum transfer cutoff 1/ro <SC Uf and the density opera- 
tors now given by p n (q) = £ Q(T J_ A {q')^(q' + q), 
and similarly for the spin-density operator p m (q) . More- 
over, the bare couplings /y in Eq. (||) should be re- 
placed by renormalized low-energy couplings gij, which 
characterize the Luttinger liquid fixed point pl|. Note 
that for ra/0 the renormalized interaction is not spin- 
rotationally invariant, so that in general g nm 7^ 0. How- 
ever, for m ^ n we expect that the generic behavior 
of correlation functions (with the possible exception of 
Xllili w ) m the spin wave regime \q\ <C q m ) can be 
correctly obtained for the special case g nm = and 

fjnn 9mm — 9 ^ ^ • 

Given the effective low-energy theory, the closed loop 
theorem jl6| guarantees that all corrections to the RPA 
for the density-density and longitudinal spin-spin corre- 
lation functions cancel for small q and to. Hence Eqs. ( |7a| ) 



and ( |7b| ) are asymptotically exact if we replace the bare 
quantities /o, Iq and Sq by the corresponding renormal- 
ized quantities g, I and 5. In particular, the existence 
of a propagating longitudinal spin mode with velocity 
v m = vfV& <C vf is a robust result, and not an artifact 
of the RPA. 

Due to the linearized energy dispersion and the irrele- 
vance of scattering processes with large momentum trans- 
fers, the single-particle Green function G a {x,r) can be 
calculated exactly using bosonization in real space and 
imaginary time. For max{|x|, Ui|r|} S> ro we obtain 



G<t(x,t) 
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Tln/2 
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[ax + iu n 7"] V 2 [aiE + IVmT] 1 / 2 



(14) 



Here rji = \(K t + K' 1 - 2), with K n = [I + I]- 1 / 2 and 
K m = [2(1— l)] -1 / 2 . Note that the anomalous dimension 
rj m of the spin channel diverges for I — > 1. This singu- 
larity is also found directly from the universal Luttinger 
liquid relation |17| Xi = ^.Kijixvi together with the above 
RPA results for Xm and v m . We note that an analogous 
scenario has recently been found for the charge channel 
of the Id t-J model in the vicinity of the phase separation 
instability Q. 

Finally, let us consider the transverse spin-spin cor- 
relation function xtli x ' T )i which, due to the linearized 
energy dispersion, can also be calculated for large x 
and t by bosonization. Following Ref. |l^ we obtain for 
max{|x|,u m |r|} > r 



-1 



2),„ 



E 



e iot{k^-k L )x 

ax + iv m T} 2 

(15) 



The leading diagrams taken into account in Eq. dl5| ) 
are shown in Fig. |^; they contain the ladder diagrams 
of Fig. ^] as a subset, but include in addition self- 
energy corrections, screening bubbles, and complicated 
vertex corrections. It is important to realize that Eq. 
( |l5| ) can only be used to obtain the Fourier transform 
Xti(q,iu) = J dxdr e~ l ( qx ~ UT " l XM( x i T ) f°r wavevectors 
close to ±(fcf — ki), i.e. for \q =p (fe-f — ki)\ < q m - In 
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FIG. 3: Dispersion of the longitudinal and the transverse spin 
excitations. The dashed line indicates that only for \q\ < 
q m we expect transverse spin waves to be well-defined. The 
triangle touching the horizontal axis at fef — fcj, is the regime 
where the bosonization result (|l^) for Sn(q, uS) can be trusted 
and yields a finite weight. The intensity of the shading is 
proportional to the magnitude of Sn(q,u>). 



the spin wave regime |g| <C q m the transverse spin-spin 
correlation function cannot be calculated using abelian 
bosonization with linearized energy dispersion, because 
(i) the ladder approximation suggests that the spin wave 
dispersion depends on the nonlinear terms of the energy 
dispersion, and (ii) the existence of spin waves follows 
from the spontaneous breaking of spin-rotational invari- 
ance, so that their dispersion cannot be obtained using a 
method which explicitly violates this symmetry. On the 
other hand, for \q =p (fcj — < q m the Fourier trans- 
form of Eq. ( |i~5| ) yields an accurate approximation for the 
transverse dynamic structure factor Sn(q,u>). For u > 
we obtain 

s h(q>v) =C m Q(uj - v m \\q\ - fc T + fcj) 



X [u> - v m (\q\ - fc T + fc x )] 
x [uj + v m (\q\ -fc T 



2)J m -l 
2r, m + l 



(16) 



with C m = [4Trv m T(2r lm )r{2 + 2T lm )]~ 1 (ro/v m ) 4 ^. The 
region where Sfi(q, w) is finite represents the Id Stoner 
continuum. The complete picture of low-energy spin ex- 
citations is depicted in Fig. || and is qualitatively quite 
similar to its 3d counterpart However, in Id there 
is no Landau damping and the structure factor shows 
anomalous scaling associated with broken spin-rotational 
symmetry of a Luttinger liquid phase. 

For an outlook from a renormalization-group perspec- 
tive, we note that while the ferromagnetic ground state 
is stabilized by nonlinear terms in the energy dispersion 
close to the Fermi points, the flow of the corresponding ir- 
relevant couplings is not accessible within the usual field- 
theoretical RG |lC| ]. However, using modern formulations 
of the RG 1 20 based on Wilson's idea of eliminating de- 
grees of freedom and rescaling, it should be possible to 
examine the subtle role played by irrelevant couplings in 
stabilizing a ferromagnetic ground state in Id. 

conclusion, we presented the effective low-energy 



of their properties only depend on the effective Stoner pa- 
rameter /, i.e., on the distance 5 = 2(1 — 1)/I -C 1 from 
the ferromagnetic instability. Neutron scattering experi- 
ments should be able to test our predictions for spin-spin 
correlation functions. Furthermore the propagating lon- 
gitudinal mode with small velocity v m cx 5 1 / 2 and large 
residue Z m oc <5 _1 / 2 dominates some thermodynamic 
quantities, for example through the divergence of the 
uniform spin susceptibility, \ m oc Z m /v m oc S^ 1 . The 
discussed features of the weakly ferromagnetic regime 
should be accessible in specially designed organic poly- 
mers Q, for which the effective Stoner parameter / can 
be controlled by adjusting the density via external gate 
voltages. Our predictions are also relevant to semicon- 
ductor quantum wires which are believed to show spon- 
taneous ferromagnetism |J 0, || . 
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